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Unified picture for the colossal thermopower compound FeSb2 
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We identify the driving mechanism of the gigantic Seebeck coefficient in FeSb 2 as the phonon-drag 
effect associated with an in-gap density of states that we demonstrate to derive from excess iron. 
We accurately model electronic and thermoelectric transport coefficients and explain the so far ill- 
understood correlation of maxima and inflection points in different response functions. Our scenario 
has far-reaching consequences for attempts to harvest the spectacular powerfactor of FeSb 2 . 
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Introduction. Thermoelectrics[l] hold great promise 
for sustainable energy solutions [2] and technological ap¬ 
plications such as Peltier coolers, heat pumps[3], micro¬ 
scopic generators [4] and probes for quality control of solid 
state devices[5]. In recent years, more and more inter- 
metallic compounds with narrow gaps have been found 
to exhibit large thermoelectrical effects. Among them, 
the iron antimonide FeSb2[6-9] is the most prominent. 
In fact, it boasts the largest thermoelectric powerfactor 
ever measured [8-11]. The thermopower assumes its max¬ 
imum at around 15 K, making FeSb 2 the prime candidate 
for thermoelectric cooling devices at cryogenic tempera¬ 
tures. The origin of the large thermoelectrical response 
is, however, unknown to date. While electronic correla¬ 
tion effects have been advocated as main benefactor [7- 
14], other works suggested the so-called phonon-drag ef¬ 
fect to cause the colossal thermopower[15, 16]. 

Besides the large magnitude of thermoelectricity, sev¬ 
eral transport quantities of FeSb 2 exhibit a pronounced 
temperature dependence[8-14]. In particular there is a 
correlation of features in different response functions[12, 
14, 17]: Extrema in the Nernst coefficient correspond to 
inflection points in both the Seebeck coefficient and the 
resistivity. Likewise, at these characteristic temperatures 
features appear in the magnetoresistance and the Hall 
coefficient. This intimate linkage between electrical and 
thermoelectrical quantities heralds a connection to the 
mechanism behind the colossal Seebeck effect. Yet, in 
all, so far no convincing comprehensive physical picture 
for this distinctive behaviour has been proposed. 

Here, we unveil the origin behind the characteristic 
temperature scales, as well as the mechanism for the 
colossal thermoelectricity in FeSb 2 . We report a min¬ 
imal model consisting of an effective one-particle elec¬ 
tronic structure with two crucial ingredients: (1) nar¬ 
row in-gap states and (2) an electron-phonon coupling 
mechanism. Our approach captures all transport phe¬ 
nomena on a quantitative level and pinpoints a new 
scenario for thermoelectricity in FeSb 2 , from which the 
phonon-drag associated with an in-gap density of states 
emerges as the main driver. Finally, using first princi¬ 
ples calculations, we identify a diminutive (intrinsic) off- 


stoichiometry, Fei+ 2 ,Sb 2 _a ;5 as likely origin of these in¬ 
gap states, in line with the large sample-dependence of 
experimental observables possibly due to different crystal 
growth techniques[6-14]. 

Theoretical framework. Employing Boltzmann’s the¬ 
ory in the relaxation time approximation[18], we study 
a variety of electrical and thermoelectrical observables. 
For the latter, we explicitly include contributions from 
the phonon-drag effect: A thermal gradient VrT pro¬ 
duces an out-of-equilibrium population in the phonon 
Brillouin zone. In case of acoustic phonons, at low tem¬ 
perature the effect is an increase in the population for 
wave vectors opposite the direction of the thermal gradi¬ 
ent. This will induce a heat flux. However, scatterings 
will drive this phonon population back towards thermal 
equilibrium. One mechanism responsible for the latter is 
the scattering with electrons, a process in which linear 
momentum is transferred to the electronic subsystem: 
the phonon-drag effect. The correct, yet cumbersome, 
way to treat these scattering events is to add a term 
to the collision integral in the Boltzmann equation. In¬ 
stead, we here propose, as an elegant short-cut, to treat 
the discrete scatterings as a continuous injection of linear 
momentum. It is easy to show that for isotropic acous¬ 
tic phonon branches, the transferred linear momentum 
is proportional to VrT. Thus an effective description of 
the phonon drag is achieved by adding a term oc VrT to 
the semi-classical equation of motion of the electrons 

k = -|[E-kf xBJ-aVrT (1) 


which states that the time derivative of the electron wave 
vector k is due to the sum of the external electrical E 
and magnetic B fields, and the contribution from the 
phonon drag {h is the reduced Planck constant, e the 
electron’s charge and r its position at time t satisfying 
f = (k), where Si is the energy dispersion of the i- 

th band). For simplicity, we shall use a band and momen¬ 
tum averaged phonon-coupling a. Then, the expressions 
for the phonon-drag driven Seebeck and Nernst coeffi¬ 
cient can be readily computed. For instance 
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FIG. 1. [Colour online] Modelled transport properties of FeSb 2 with one in-gap state compared to experiment. [14] (a) density 
of states (left) and chemical potential p as a function of temperature (right), (b) electrical conductivity a on a logarithmic 
scale, consisting of valence (V.B.), conduction (C.B.) and in-gap band (B.l) contributions, (c) Seebeck coefficient, (d) mag¬ 
netoresistance (MR), (e) Hall coefficient Rh, and (f) Nernst coefficient. The (gray) vertical line indicates the temperature at 
which the chemical potential passes the in-gap state i?i. 


where / is the Fermi function, the relaxation time, Ni 
the band degeneracy, and a the total conductivity. 

Evaluating the transport functions requires the knowl¬ 
edge of the dispersions £i. We choose parabolic conduc¬ 
tion and valence bands with equal mass. However, the 
latter are not enough for a complete description. Com¬ 
paring the experimental resistivity of FeSb 2 to the ac¬ 
tivated behaviour of a perfect semiconductor (see Sup¬ 
plemental Material), we notice that the shoulder around 
lOK is a clear fingerprint of defect states inside the gap 
(as also reported in Refs. 8, 10, and 19). We therefore 
explicitly include in-gap states in our treatment and con¬ 
sider, as a practical approximation, pairs of spherically 
symmetric dispersions: (fc) = Ei ± Wi/ 2 cos{'^^/ 2 k), 

where the reciprocal lattice vector k lies in a Brillouin 
zone of radius k and Ei is the energy position and Wi 
the width of the band. 

In the following, it will be important to distinguish be¬ 
tween the temperature dependence and the amplitude of 
the response functions. It can be proven[18] that each 
transport properties can be expressed by an amplitude- 
factor and a function rj that contains all non-trivial tem¬ 
perature dependence. For example, the contribution to 
the phonon-drag thermopower of an in-gap state can be 
written as: 

^ aWfhenN, ^ (E,-pi{T) W, \ 

192h^akBT ksT ’ 2kBTJ ' 

Similar expressions hold for all transport properties, with 
different amplitude prefactors (depending on Wi, Ni, Ti, 
ki) and characteristic rj-functions that all depend on the 
same two variables: and ry-functions 

can be classified into two categories according to their 
symmetry with respect to their first argument. For exam¬ 
ple, the conductivity, which does not distinguish between 
electrons and holes will have the same amplitude whether 
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TABLE I. Schematic representation of the even/odd character 
of the ry-functions for in-gap state derived transport. 


the chemical potential is above or below the band cen¬ 
tre, hence rja is even. Conversely, the Hall coefficient 
has to change sign when the chemical potential is at the 
band centre; vibh is odd in its first argument. The cate¬ 
gory of an ? 7 -functions thus reflects the sensitivity of the 
transport property to the type of carriers. A full de¬ 
scription can be found in the Supplemental Material[18], 
while a quick reference is reported in Table. I. It is re¬ 
markable that while the electronic Seebeck coefficient is 
sensitive to the carrier type, the phonon-drag driven one 
is not. As a consequence the relative symmetry between 
S, N, Rh and a can be used to distinguish purely elec¬ 
tronic from phonon-drag driven thermoelectricity. Inter¬ 
estingly, we also find that all yy-functions are only slowly 
varying with their second argument. In particular the 
qualitative shape and the symmetry are inert, and even 
the peak widths are only weakly dependent on 2 ^^- 
Consequently the temperature profile of all the transport 
properties deriving from in-gap states is entirely deter¬ 
mined by the position Ei of the band with respect to the 
chemical potential and the absolute temperature. All 
other parameters merely rescale the amplitude in a T- 
independent fashion. 

Temperature dependence of transport properties. We 
now allow for one state £i{k) inside the charge gap, 
A = SOmeV, to account for the shoulder in the conduc¬ 
tivity (we find that Ei = 24.3meV yields good agree¬ 
ment, see Fig. 1(a)). All other parameters {Wi, Ni, • • •) 
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FIG. 2. [Colour online] Same as Fig. 1 but for three in-gap states: Ei = 11.4, 19.8, 24.3 meV. The (gray) vertical lines are guides 
to the eye indicating maxima in the MR. For details of the fitting procedure see the text and the Supplemental Material[18]. 


merely describe the relative amplitudes for the conduc¬ 
tion/valence and in-gap state contributions, and will be 
discussed later. Here, we first focus on the temperature 
dependence of the thermoelectric response to elucidate 
the ill-understood correlations between features in dif¬ 
ferent transport quantities of FeSb 2 . We consider three 
scenarios: (1) a purely electronic description of thermo¬ 
electricity (with contributions from the VB, CB and in¬ 
gap states), (2) a dominant phonon-drag effect on va¬ 
lence and conduction carriers, and (3) a preponderant 
phonon-drag effect that couples to electronic in-gap ex¬ 
citations. Scenario (1) is known to be insufficient to ac¬ 
count for the amplitude of thermoelectricity in FeSb2[15], 
but it also yields the wrong temperature profile. Scenar¬ 
ios (1) and (2) are dominated by valence and conduction 
contributions so that the Seebeck and the Nernst coef¬ 
ficient have a very comparable temperature dependence 
and peak positions -at variance with experiments. For 
details of the hence discarded scenarios see the Supple¬ 
mental Material[18]. 

This leaves us with scenario (3), a preponderant 
phonon-drag effect linked to in-gap states, that we ad¬ 
vocate in the following. In this picture, the shapes of all 
transport functions are uniquely determined by a single 
parameter: the position Ei of the in-gap state (see the 
discussion of Eq. (3)). As can be seen in Fig. 1, this 
setup correctly predicts the characteristic peak temper¬ 
ature, as well as the line shapes of the major feature in 
all five transport functions. 

The good agreement notwithstanding, there are devi¬ 
ations above 15K. To understand this issue, we scruti¬ 
nize the experimental magnetoresistance (MR). We no¬ 
tice that FeSb 2 displays a strong response to a magnetic 
field only in the temperature range where our in-gap band 
gives a contribution to the conductivity, see Fig. 1. We do 
not provide an expression for the unconventional linear- 
in-B MR of this material, as the underlying mechanism 
is unresolved. Nonetheless, we can argue that the MR 
cannot distinguish between electron and hole transport 


and has an even-type ? 7 -function. Consequently, the MR 
peaks at the same temperature as the conductivity (the 
peak width might be different). In Fig. Id we there¬ 
fore compare the MR to the in-gap-state conductivity 
and immediately notice that one in-gap band gives only 
one peak in the MR. However, in the experimental data 
there are in total three features. In Fig. 2, we there¬ 
fore refine our model to include three in-gap states (note 
that a single state with an energy differentiation of its 
density of states or group velocities can play the same 
role). All peaks and their correlation in different trans¬ 
port functions are then captured with only three param¬ 
eters: the energy positions Ei of the in-gap states. To 
account for relative amplitudes (see below), our model in¬ 
cludes two further parameters that we extract from the 
ratios of the Nernst coefficient peak amplitudes. Then, 
all five transport coefficients follow accurately: In Fig. 2 
the peak positions and line widths of all features are 
captured [20]. This pinpoints the preponderant influence 
of in-gap states and their role in the phonon-drag picture 
as the main driver for the physics of FeSb 2 . 

Amplitude of transport coefficients. While we focused 
so far on the temperature profile of transport coefficients, 
the modelling in Fig. 2 also neatly accounts for the in¬ 
dividual amplitudes. We now detail the parameters that 
quantitatively determine the response functions. The 
phonon-drag constant can be extracted completely from 
experiments as a » 3 • l{f^^s~^K~^ = 229.2 fcs/Ti at 
T = 10K[I8]. 

Since the temperature dependence of the transport 
functions is uniquely controlled by the positions of the 
in-gap states, every transport coefficient only contributes 
one independent value for the determination of the band- 
parameters, Wi, Ni, Ti. We are thus dealing with a 
largely under-determined set of equations. Therefore, 
we fix some of the parameters of our model to physi¬ 
cally reasonable values and then check all remaining ad¬ 
justable parameters for consistency: We know that W 
has to be of the order of the unity, so we assign Ni = 1 





































4 


for all the in-gap bands. Then, the width of the Bril- 
louin zone for in-gap states (we assume = k) can 

be extracted from the ratio of the peak amplitudes of the 
Nernst and Seebeck coefficient, yielding k « 10“^A“^. 
Comparing k to the size of the Brillouin zone of per¬ 
fect FeSb 2 , we deduce an effective defect concentration of 
9 • 10“®/unit-cell. Next, we set = tcb = tvb = O.Sps. 
This fully determines the in-gap bandwidths for which 
we find Wi = 3meV, W 2 = 4.1 and W 3 = 5.1 meV [21]. 

With all parameters thus fixed, we now check if the 
value of a is realistic. To this end we compare the lin¬ 
ear momentum that is transferred from the phonons to 
the electrons with that dissipated by lattice thermal con¬ 
duction. We find the momentum transfer to the elec¬ 
trons -that is needed to justify the large S- to be only 
~ 0.5% of that dissipated by phonons (for details see Sup¬ 
plemental Material[18]), which is very realistic. Indeed 
the experimental thermal conductivity is characterised 
by phonon-crystal boundary scatterings [8], supporting 
such a small momentum transfer. Consequently our re¬ 
sults show that a large phonon-drag thermopower need 
not produce temperature characteristics in the thermal 
conductivity. That a minute transfer of phonon momen¬ 
tum has such a big impact on charge carriers is explained 
by the large ratio of masses of ions and electrons. 

In the Nernst coefficient, the proposed phonon- 
coupling to the density of in-gap states circumvents Sond- 
heimer cancellation[22], without requiring multiband ef¬ 
fects. Previously, it had been proposed that a strongly 
energy-dependent scattering rate r(e) causes the colos¬ 
sal Nernst signal[14]. Indeed we can assume an e/- 
fective, purely electronic picture with a scattering rate 
T(e) oc e’’, where the exponent is determined by r = 
— e/fcs X N/{Rh<j)- This yields, in our phonon-drag for¬ 
malism, a large exponent r = —(h/kB)a = —229.2, in 
reasonable agreement with Ref. 14. 

Origin of in-gap states. We now investigate the pos¬ 
sible origin of the in-gap states. The number of states 
has to be low compared to the weight of the bottom 
of the conduction band, or the chemical potential will 
not be able to cross them. Moreover the in-gap band¬ 
width must be small compared to the intrinsic charge 
gap. This suggests the states to stem from local de¬ 
fects. The strong reaction to magnetic fields (MR) then 
points towards Fe impurities. Using density functional 
theory calculations, we have modelled the behavior of 
FeSb 2 with different defects and impurities. We find 
that substituting Fe on Sb-sites, i.e. Fei+a;Sb 2 -£c, agrees 
best with the picture proposed above. Our estimate for 
the minimal defect concentration (see above) translates 
to a: « 4.5 • 10“^ [18]. Simulating the effect of iron de¬ 
fects, we show in Fig. 3 the paramagnetic band struc¬ 
ture of Fe 25 Sb 47 (i.e. x « 0.04) [23]: Indeed two impu¬ 
rity bands appear close to the conduction states. More¬ 
over, the doping coming from the defect pins the chem¬ 
ical potential (at T = 0) to the bottom of the conduc- 
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FIG. 3. [Colour online] Band-structure of Fe 25 Sb 47 , mimick¬ 
ing Fei+a;Sb 2 - 2 :. Hues of green indicate orbital admixtures 
originating from the extra iron atom. Inset: density of states 
of pure (black) and Fe-rich (green) FeSb 2 . 

tion band, in extremely good agreement with the sce¬ 
nario proposed above. This finding is supported by the 
binary Fe-Sb phase diagram[24] that may suggest Fe- 
richness, and a recent experiment demonstrating that 
even the smallest iron concentrations in Rui_ 2 ;Fea;Sb 2 
causes the appearance of the low temperature character¬ 
istics of FeSb2[13]. Further support comes from recent 
magnetic resonance experiments [25] that advocate local¬ 
ized S' = 1/2 in-gap states slightly below the conduction 
bands. The presence of Fe defects with localized states 
entails a Curie-like contribution to the magnetic suscep¬ 
tibility, which can be estimated as follows: Curie’s law 
reads y = xp,oNAfJ,1ff /with the vacuum perme¬ 
ability fiQ, Avogadro’s constant Na, an effective local mo¬ 
ment fXeff = gs\/S{S l)/r_B, with gyromagnetic factor 
gs = 2 and Bohr magneton piB- Using the experimentally 
suggested S = 1/2 spin-state[25] and the minimal defect 
concentration x yields x • T = 5 • 10“®emuK/mol. This 
is consistent with experiment, witnessing a low tempera¬ 
ture upturn in the susceptibility for most samples, with a 
magnitude smaller than 5 • 10“^emu/mol[6, 7, 9, 10, 26- 
28] for all accessed temperatures. 

Summary. Purely electronic transport properties of 
FeSb 2 , like the conductivity and the Hall coefficient, 
can be described with impressive accuracy from an ef¬ 
fective one-particle electronic structure [29] when taking 
into account the influence of defect-derived in-gap states. 
Using ab initio band-theory, we identified surplus iron, 
Fei+a;Sb 2 -a;, as the likely origin of these states. Our 
work provides a unified picture: there is only one inde¬ 
pendent transport property, from which the features of 
all five other response functions are determined in a one- 
to-one correspondence. Finally, we have pin-pointed the 
phonon-drag associated with the density of in-gap states 
as the microscopic origin of the hitherto elusive colos- 
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sal thermoelectricity in FeSb 2 . Our findings have impor¬ 
tant consequences for the quest of reducing the thermal 
conduction. Selective phonon engineering, e.g. by nano- 
structuring[3, 30], needs to take care that the phonon 
drag mechanism survives. Our scenario might also be 
relevant for understanding other narrow-gap semiconduc¬ 
tors, such as CrSb2[31] or FeGa3[32, 33]. 
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